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Abstract. We consider the numerical approximation of the spectrum of a second- 
order elliptic eigenvalue problem by the hybridizable discontinuous Galerkin (HDG) 
method. We show for problems with smooth cigenfunctions that the approximate 
eigenvalues and eigenfunctions converge at the rate 2k + 1 and k + 1, respectively. 
Here k is the degree of the polynomials used to approximate the solution, its flux, 
and the numerical traces. Our numerical studies show that a Rayleigh quotient-like 
formula applied to certain locally postprocessed approximations can yield eigenvalues 
that converge faster at the rate 2k + 2. We also derive and study a condensed nonlinear 
eigenproblem for the numerical traces obtained by eliminating all the other variables. 



We study the HDG (hybridized discontinuous Galerkin) approximation to the follow- 
ing eigenproblem: Find eigenvalues A in R and corresponding nontrivial eigenfunctions 
u satisfying 



Several HDG discretizations were introduced in [8] to discretize corresponding source 
problems. The purpose of this paper is to study the application of one such method to 
eigenproblems. The particular HDG method considered here (referred to as the H-LDG 
method in [8], but simply as the 'HDG method' in this paper) is chosen for our study 
because we now have a fairly complete theoretical understanding of its application to 
the source problem. 

Two well-known advantages of the HDG method, when applied to source problems, 
are reduced system size and flexible stabilization. The latter arises due to a transparent 
stabilization mechanism in the definition of numerical flux. It allows one to use, for 
example, polynomials of the same degree k to approximate the solution u and com- 
ponents of its flux q = —aVu. While use of these spaces would have resulted in an 
unstable mixed method, the resulting HDG method is stable, and provides optimal or- 
der approximations for all variables. To discuss the former advantage of reduced system 
size, recall the process of static condensation, which, for source problems, removes all 
interior variables to yield a "condensed" system for inter-element variables. The HDG 
condensed system, when compared to condensed systems from other DG methods, is at- 
tractive because of its smaller size and favorable sparsity patterns. One of the questions 
we address in this paper is whether such condensed systems are useful in eigenproblems. 
Note that the condensation process, by reducing the system size, also reduces the size of 
the spectrum. How much of the spectrum can be approximated despite this reduction 
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is a question answered in Section [5j There, we also derive the nonlinear equation that 
needs to be solved in the reduced dimensions to compute the spectral approximations. 

Apart from these results on the condensed eigenproblem, we prove convergence results 
for the HDG eigenproblem in Sections [3] and IH We show that the approximate eigenval- 
ues obtained by the HDG method exhibit no spectral pollution. They converge to the 
exact eigenvalues at the rate 0(h 2k+l ), under suitable regularity assumptions, if we use 
polynomials of degree at most k for all the HDG variables. We also show, under similar 
assumptions, that the gap between the corresponding discrete eigenspace and the exact 
eigenspace in L 2 (Q) converges at 0(h k+1 ). Roughly speaking, this shows that the rate 
of convergence of eigenf unctions in L 2 (Q) is optimal for the HDG method. 

These results compare favorably with many other DG eigenvalue approximations [21 
[T3"] . The unified presentation of [2] facilitates comparison. They show that many tra- 
ditional Hermitian DG methods approximate eigenvalues at the rate 0(h 2h ). For non- 
Hermitian DG methods, they find that the eigenvalue convergence rate is even poorer, 
in general no better than 0(h k ). The HDG method, which can be considered to fall in 
the Hermitian class, yields eigenvalues that converge faster, when compared to both the 
Hermitian and the non- Hermitian DG methods considered in [2]. However, let us note 
that the convergence rate of HDG eigenvalues (or, for that matter, any DG eigenvalues) 
do not compare favorably with the 0(h 2h+2 ) convergence rate of the mixed hybridized 
Raviart-Thomas (HRT) method [SI El US]- Our analysis also points to other differences 
when comparing the HDG and HRT eigenproblems. For example, comparing Theo- 
rem [53] below with [91 Theorem 3.2], we note that the extent of the spectrum recovered 
by the condensed system may be shorter (up to 0(1/ h)) for the HDG method in com- 
parison to the HRT method (which is up to 0(l/h 2 )). 

The method of convergence analysis in this paper is motivated by the many early works 
that developed abstract approaches to analyze approximation of eigenproblems [6l [TTl[T6"l 
[XT] , and in particular, the application of the abstract theory to DG methods in [2]. The 
critical new tool that helps push the analysis forward in the HDG case, and yield better 
convergence rates than [2J, is the projection operator of [TU]. The projection allows our 
analysis of HDG eigenvalue errors to proceed along the lines of similar analyses for mixed 
methods [5]. A few important differences arise due to the fact that the HDG projection 
possesses only a weak analogue of a well-known commutativity property. The analysis 
in the second aspect of this study, involving the condensed system, is motivated by our 
previous such analysis [9] for the HRT method. 

Reports of extensive numerical experiments are provided in Section [61 Of particular 
interest is a local and inexpensive postprocessing presented there. The postprocessed 
eigenvalues seem to converge at 0(/z 2fc+2 )-rate thus making the HDG method competitive 
with the mixed method. The intuition behind the construction of this postprocessing is 
inherited from our previous experience with the mixed method [9], where we proved that 
the postprocessed eigenfunctions have better convergence rates. However, the argument 
used there to provide a rigorous proof does not seem to extend to the HDG method. 

In the next section, we introduce the HDG eigenproblem and essential notations used 
throughout the paper. In Section [31 we show that there is no pollution of the spectrum 
when approximating it by the HDG method. This result is improved in Section HI where 
we establish (in Theorem 14. ip the convergence rate of HDG eigenvalues. Section \5\ 
investigates the condensation (or hybridization) of the HDG eigenproblem to obtain a 
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smaller condensed nonlinear eigenproblem for the interface variable otherwise known as 
the numerical trace. The condensed nonlinear eigenproblem is given in Theorem 15.31 
and a closely related linear eigenproblem is investigated in Theorem 15.41 We conclude 
in Section [6] with the results of our numerical studies. 

2. The HDG source and eigenvalue problems 

Consider the Dirichlet boundary value problem (rewritten as a first order system) of 
finding qf £ H(div, ft) and vJ £ L 2 (ft), given any "source" / in L 2 (ft), such that 

q f + a Vu f = on ft, (2a) 

V • q f = f on ft, (2b) 

u f = on <9ft. (2c) 

All functions, unless explicitly stated otherwise, are real-valued in this paper. Through- 
out, ft C R n is a polyhedral domain (n > 2), a : ft t— > M. nxn denotes a variable ma- 
trix valued coefficient, which we assume to be symmetric and positive definite at all 
points in ft. To facilitate our analysis, we introduce notation for the "solution operator" 
T : L 2 (ft) i — y L 2 (ft), which is defined simply by 

Tf = uK (3) 

It is well known that T is compact and self-adjoint. Its spectrum, denoted by cr(T), 
consists of isolated points on the positive real line accumulating at zero. Clearly, there 
is a one to one correspondence between the eigenvalues of (jrj and those of T. Indeed, 
fj, is an eigenvalue of T if and only if n = 1/A for some A satisfying ([I]) . 

2.1. The source problem. The HDG method provides an approximation Th to T. To 
understand this approximation, we first describe the HDG source problem and introduce 
known results we shall use later. Afterward, we will present the HDG eigenvalue problem. 

The HDG method defines a scalar approximation Uh to u and a vector approximation 
cfh to q in the following spaces, respectively: 

Wh = {w : for every mesh element K,w\k £ Pk{K)}. (4) 

Vh = {v : for every mesh element K,v\k £ Pk(K) n } (5) 

Note that functions in these spaces need not be continuous across element interfaces. 
Above and elsewhere, we use Pk{D) to denote the space of polynomials of degree at 
most k > on some domain D. The subscript h denotes the mesh size defined as the 
maximum of the diameters of all mesh elements. 

For any (scalar or vector) function q in Vh or Wh, the trace q\p is, in general, a 
double-valued function on any interior mesh face F = dK + D dK~ shared by the mesh 
elements K + and K~ . Its two branches, denoted by [q]x+ and are defined by 

[q]K±(x) = lim^o q(x — e[n]K±) for all x in F. Here and elsewhere, n denotes the 
double-valued function of unit normals on the element interfaces: on a face F C dK, 
its branch [n\K equals the unit normal on OK pointing outward from K. For functions 
u and v in L 2 (D), we write (u, v)p = j D uv dx whenever D is a domain of M n , and 
(u,v)d = J'd uv dx whenever D is a domain of M n_1 . To simplify the notation, define 

{v,w)sr h = ^ ( v i w )k and (v,w) 9 ^ h = ^ ( v i w )dK, 
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where in the latter, we understand that for double valued v and w, the integral (v, w)qk 
is computed using the branches [v]k and [w]k from K. For vector functions v and w, 
the notations are similarly defined with the integrand being the dot product v ■ w. 

In addition to the spaces Vh and Wh introduced above, our method also uses one other 
discrete space Mh, consisting of functions defined on the mesh faces (or mesh edges if 
n = 2), namely 

Mh = {/i : for every mesh face F, [i\p G Pk(F), and if F C dQ, n\p = 0}. (6) 

The HDG method defines the approximate solution ul, the approximate flux q^, and 
the numerical trace r)l, as the functions in Wh, Vh, and M h , respectively, satisfying 

0, for all re V h , (7a) 

(f,w),? h for all weW h , (7b) 
for all fj, G M h , (7c) 

where c = a~ x and q\ is a double- valued vector function on mesh interfaces defined by 

fh = € + T ( u l-V f h )fi. (8) 

Note that this defines all branches, i.e., on the boundary OK of every mesh element K, 
the value of the branch of q[ from K is [qf^K = [%\ K + [ r W ~~ Vh) [™\k- Here r is 

a non-negative penalty function, which is also double-valued on the element interfaces 
and [t]k above denotes the branch of r- values from K. For simplicity, we assume that 
any branch of r is a constant function on each mesh edge. It is proved in [S] that the 
system (J7]) is uniquely solvable if tk is positive on at least one face of K for every 
element K. This unique solvability result is assumed throughout this paper. Given any 
/ in L 2 (Q), the component w{ of the unique solution of (JTJ) is used to define the discrete 
version of the operator T in namely 

T h f = u{. (9) 

We will need a projection IIh(q, u), into the product space Vh x Wh, originally designed 
in [10]. Its domain is a subspace of H(div, Q) x L 2 (Q) consisting of sufficiently regular 
functions, e.g., H(dxv, Q) fl if s (fi) n x H s (il) for s > 1/2. When its components need 
to be identified, we also write n h (q,u) as (n%q, II^u) where n^q and n^u are the 
components of the projection in Vh and Wh, respectively. (Despite this notation, note 
that JJ^q depends not just on q, but rather on both q and u. The same applies for 
II^u.) The components are defined by 

(q,?)K for all r eP k _ x {K) n , (10a) 

(u,w) K for all w G P k -i(K), (10b) 

(q-n + ru,fi)F for all /i G Pfc(F), (10c) 

for all faces F of the simplex K. Let s u ,s q G (1/2, k + 1]. We recall the following 
approximation property, proved in [TD], Theorem 2.1] for integer values of s u ,s q , and 
extended to remaining values of s u , s q G (1/2, k + 1] in [7]: 



-(q[, Vw) Sh + (f h ■ w )a? h = 
{^,ql-n)a^ h = 



(n h q,r) K = 
(n h v u,w) K = 
{n h / q-n + TlI h ¥ u,n) F = 
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For all (g, u) G if(div, fi) n H s "(VL) n x # s «(fi), 

||i7£g-g|U < C/# Igl^jo + C/i^ri |«| H ^ (J0 (Ha) 

\\n%u-u\\ K < Ch% \u\ H su {K ) + C ^\q\ H s q{K) . (lib) 



Here = maxr\gK\F max , where F max is a face of K at which t\qk is maximum. And 
is the diameter of K. We also denote h = max^ hx- Above and throughout we use C to 
denote a generic constant independent of the mesh element sizes and the stabilization 
parameter r. The following error estimate is known [7J [10] . 



Theorem 2.1 (see [7]). Let the exact solution (q f , u f ) of © be in #(div, fi)nif s (fi) n x 
H S (Q) for some s > 1/2. T/ien, 

ll«-«£lk <cil«-^r«lk + &rC||g-ii^lc, (12) 

\\q-qt\\ c <2\\q-nlq\\ c , (13) 

1 /2 

where b T = max{l + HrT^ + hx/T^ ax : il" G and ||g || c = (eg, q)g^ h - 

2.2. The eigenproblem. The HDG discretization of the eigenproblem ([1]) defines an 
approximation to the eigenfunction Uh G Wh, an approximation the eigenflux q\ G V^, 
and an approxmation to the eigenfunction trace rjh G M^, as a nontrivial set of functions 
satisfying 

{cq h ,r)& h - (u h , V-r)sr h + {Vh,r -n)a^ = 0, for all rG 14, (14a) 

-(&, Vto)^ + (g h • n, iw)a% = K (u h , w) :3 r h for all w E W h , (14b) 
(/x, g ft • n) a s h = for all /j G M h , (14c) 

where g^ is defined by % = qh + T\Uh — rjhjn, cf. (jSJ). Here, A^ G K. is the corresponding 
discrete eigenvalue. 

The unique solvability of the source problem ([7]) implies that Xh is nonzero. One can 
easily verify that \ih is an eigenvalue of Th if and only if ^ — 1/ A^ for some A^ solving 
(I14p . Moreover, A^ is positive as can be concluded from the next lemma. 



Lemma 2.2. 7) t is self-adjoint and positive definite in L 2 (Q). 

Proof. To show that (f,T h g) = (T h f,g), where (•, •) denotes the L 2 (f2)-inner product, 
we calculate as follows: 

(/, T h g) = -(g/, V<)^„ + (q[ ■ ft, <)a^, by (|7b]) with w = T h g = u 9 h 

= (V • ql , ul), h + ((gt - <f/) • n, <U by integration by parts 

= (cglf, + g/ ■ + - €) ■ * by © with f = g*/ 

= g/),% - (g/ • - ^>ay h + (g{ • rf, < - rfymr,, by (JTcJ) with [i = rf h 

= (cgT, gf)^ + (tK - « - Vl))a? h by ©. 

The last expression is symmetric in / and g and is non- negative if / = g. This proves 
that Th is self-adjoint and positive semidefinite in L 2 (Q). As already noted previously, 
zero is not an eigenvalue of Th, hence Th is positive definite. □ 
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3. Convergence of the spectrum 

The convergence of the discrete eigenvalues to the exact ones is proved by establishing 
convergence of T h to T in operator norm. (Recall that T and T h are defined in ([3D and 
respectively.) Such operator convergence was used as the basis for the early analyses 
of spectral approximations using conforming methods [61 [161 ELT] • It h as a l so been used 
to analyze approximations of eigenvalue problems using older discontinuous Galerkin 
methods (like the interior penalty method) [2]. To apply this technique to the HDG 
eigenproblem, we need the following basic result. 

Theorem 3.1 (Operator convergence). Suppose there is an s > 1/2 such that any 
solution (qf,u-f) of (T5[) satisfies 

\W f \\H'(a) + H^llff-cn) < C||/IU2(n), (15) 

for all f E L 2 (n). Then 

\\T-T h \\ < c T /i min(s ' fc+1) (16) 
where || ■ || denotes the L 2 {Vt)- operator norm and c T = C max{l + h 2 K (r^)' i + (r]£ ax )~ 2 : 

K e s? h y/ 2 . 

Proof. The convergence results for the HDG source problem imply 

\\Tf - T h f\\ L 2 {n) < C\\u - II™u\\ + b T C\\q- W h q\\ c , by TheoremEU eq. (JED, 
< Cc T h mia(s > k+1 \\q\ HS{n) + \u\ H . m ), by (P). 
Hence the result follows from (fl5l) . □ 

Note that assumption f[T5|) is a regularity assumption that holds, for example, when 
s = 1 and Q is a polygonal domain (with no cracks) in IR 2 . 

By virtue of Theorem 13. II and the well-known consequences of operator convergence [TJ 
[12l ITT] , we conclude that the spectrum of approximates that of T, i.e., there is no 
"pollution" of the spectrum when it is approximated by the HDG method. To formulate 
the statement of this approximation precisely in a form we can use later, let us recall 
some standard terminology. The "gap" between two subspaces X and Y of L 2 (Q) is 
defined by 

dist(x,F) dist(y,X) 

d(X,Y) = sup — = sup — . 

x&x \\x\\ yeY \\y\\ 

Here and throughout || ■ || denotes the L 2 (f2)-norm. 

Now, suppose fj, is an eigenvalue of T of multiplicity m and let r be a positively 
oriented circle, contained in the resolvent set of T, centered at /i, and enclosing no other 
eigenvalue of T. Define two operators, E£ and E r , both on L 2 (Q), by the following 
integrals over r in the complex plane: 

E r h =±-l{z- T^dz, E r = ±<f(z- TY l dz. 

Hereon, we omit the superscript r in E? and E r as r will always be taken as stated 
above. A well-known result is that E is a projection onto the eigenspace of T correspond- 
ing to the eigenvalue fi. The discrete analogue of this result for Eh appears in results 
collected below. The collection summarizes a few consequences of operator convergence. 
The arguments proving these consequences are standard [2l[T6j[TT], and since they apply 
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to in the HDG context with few modifications, we shall not repeat them. We use R(A) 
to denote the range of any operator A and C T to denote generic constants independent 
of h, but dependent on r. 

Corollary 3.2 (Convergence of eigenvalues and eigenfunctions) . Let fi and T be as 

above and let the assumption of Theorem \3. II hold for some s > 1/2. Then, there exists 
hp > such that for all h < hp the following statements hold: 

(1) There are exactly m eigenvalues of T^ within F , which we count according to 
multiplicity and enumerate as fih,i, H>h,2, ■ ■ ■ > ^h,m- 

(2) The operator Eh is a projection onto the span of the eigenfunctions ofTh corre- 
sponding to all the eigenvalues fi^j for j = 1, . . . , m. 

(3) The operator E^ converges to E as h — > and 

\\E - E h \\ < C\\T - T h \\ < C T h min{s ' k+1) . (17) 

(4) The exact and discrete eigenspaces (of fi and {aa/ij}^L 1; respectively) are R{E) 
and R(Eh), respectively. The discrete eigenspaces converge in the sense that 

5{R{E),R{E h )) < C T /i min(s ' fc+1) . (18) 

(5) If in addition, the eigenfunctions of \i have a higher regularity index, i.e., if 

H^lk^n) + ll^H^cn) < C\\f\\ L 2 m , \Jf g R(E), 

with s M > s, then (1171) can be refined to 

\\(E - E h )\ R[E) \\ < C\\(T-T h )\ R{E) \\ < C^ min(s - fe+1) . (19) 
and consequently (JT8J) can be revised to 5(R(E), R(E h )) < C r /i min(s, " fc+1) . 

4. Rate of convergence of eigenvalues 

In this section, we prove that under favorable regularity conditions, the HDG eigen- 
values converge at the rate 0(h 2k+1 ) when we use polynomials of degree at most k > 
for all variables. To do so, we begin with the setting of Corollary 13.21 and refine a few 
estimates through a duality argument. Accordingly, we keep the same notations as in 
Corollary I3.2[ and tacitly assume throughout this section that the assumptions in the 
corollary hold. In particular, recall that R(E) is the eigenspace of T corresponding to 
an eigenvalue /i and fihj are the discrete eigenvalues near /i. 

Theorem 4.1. Suppose there is an > 1/2 such that 

\\q f \\H^(n) + ll^ll^+i^) < C\\f\\ L 2 {n) , Mf G R{E). (20) 

Then there is an /i M > such that for all h <h il , 

\p ~ HhJ < Ch m ' m ^' k+1)+m ' m ^' k \ (21) 

This is the main result of this section, and in the remainder of this section, we prove it. 
As we shall see, we are able to apply the needed duality techniques for this proof, thanks 
to properties of the projection operator LT h (q, u) of [TO], recalled in Subsection 12.11 Given 
any eigenfunction e G R{E), let us denote its corresponding flux by 

q e = -«Ve. (22) 
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Then H™e is the component of 77/i((f e ,e) in Wh- Define 

J h e = E h n™e 

for all e in R(E). With these notations, we begin the analysis with the following two 
lemmas. 

Lemma 4.2. Eh is a self-adjoint operator in L 2 (Q). 

Proof. By Lemma 12721 is self-adjoint. It is well-known that the spectral projection of 
any normal operator is self-adjoint [15]. Since Eh is the spectral projection of Th, the 
lemma follows. □ 

Lemma 4.3. Suppose (120]) holds. Then, there exists an ho > such that for all h < ho, 
the operators Jh : R(E) i — > R(Eh) and Eh\ R ^ ■ R(E) i — > R(Eh) are bijections, and 
there are h-independent constants C T ,j such that 

a,i||e|| < ||4e|| <C r ,2||e||, (23) 
a, 3 ||e|| < \\E h e\\ < \\e\\, (24) 

for all e in R(E). 

Proof. By item H] of Corollary 13.21 the gap between R(E) and R{Eh) becomes less than 
one when h is small enough, and consequently dim(R(E)) = dim(R(E h )) (see, e.g., [HJ 
Lemma 221]). Therefore, to prove the stated bijectivity, we only need prove injectivity. 

In preparation, we recall that by Lemma |4~2"| Eh is self-adjoint, and by Corollary l3.2[|2"j ). 
Eh is a projector. Hence Eh is an orthogonal projector. Orthogonal projectors have unit 
norm, hence 

\\E h v\\<\\vl WveL 2 (Q). (25) 
A second preparatory inequality we need is 

|| e - i7f e|| < C T h r \\e\\, Ve G R(E), (26) 

with r = min(s /J , k + 1). This is a consequence of (lllb|) . by which 

II e - 11% e\\ < C T h r (\e\ H s^n) + \q e \ H ^(a)), 

where q e is as in ( 122|) . Thus, ( )26|) follows from the regularity assumption (|20|) . 
Let us now prove (123|) . beginning with the lower bound. 

\\J h e\\ > ||e||-||e-^i7^ e || 

= ||e|| - \\(Ee - E h e) + E h (e - It^e)\\ since e G 

> ||e|| - \\{E-E h )e\\ - \\e-TI%e\\ by (T25]) 

> (1 - C T h r )\\e\\ by (ESJ) and JTS}. 

Therefore, the lower bound follows by choosing small enough h. The injectivity of Jh 
is an obvious consequence of this lower bound. The upper bound of (123 p immediately 
follows by combining (I26p . 

Wh v e\\ < (1 + C T h r )\\e\\ Ve G R(E), 

with (J2SD- 

To finish the proof, note that the upper bound in is already proved in ( l2"5"j) . The 
lower bound in ( 1!M|) (and the consequent injectivity of Eh on R(E)) follows by an simpler 
argument similar to the above. □ 
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Proof of Theorem J^.l. As a first step, we define two finite dimensional operators whose 
eigenvalues are \i and fihj- Let h be so small that we can conclude by Lemma [4.31 that 
J^ 1 : R(E h ) H- R(E) exists. Set T = T\ R r E ) and T h = J^ l T h J h \ R ^ E y Both the operators 

f : R(E) M- R(E) and f h : R(E) H- R(E) 

are finite dimensional. T has /i as its (only) eigenvalue of of multiplicity m. Moreover, it 
is easy to see that are the eigenvalues of T h . Hence, by the Bauer-Fike theorem [3J, 

\n-H h j\<C\\f-f h \\. (27) 

The remainder of the proof bounds the above right hand side appropriately. 
To this end, let / 6 R(E), and consider (T — T h )f. Then, 

C T>l \\{f-f h )f\\ < \\J h {f-f h )f\\ by d2D of LemmaSJ 

= \\E h n%Tf-T h E h n™f\\ 

= \\E h {n?Tf - T h n™f) || as T h and £ h commute. (28) 

We bound the norm in ([281 by duality, as follows. By Lemma \A. 3 1 any in R(Eh) can 
be written as E^g for some g in R(E). Therefore, 

\\E h (IIZTf-T h nZf)\\= sup ( E >r( n kTf-T h n^f),g h ) 

9h&R{E h ) \\9h\\ 

= gup {E h (II?Tf-T h II?f),E h g) 
geR(E) \\ E h9\\ 

< — — sup v fe fc - (29) 

^t,3 g&R{E) \\9\\ 

Note that to obtain the numerator above, we used the self-adjointness of Eh given by 
Lemma 14. 2[ while to obtain the denominator, we used Lemma 14.31 

It will now be convenient to split the numerator in (f2"9"j) into several terms. With / 
and g in R(E), we write 

(n^Tf - T h II™f, E h g) = (II^Tf - T h f, E h g) + (T h (f - Il^f), E h g) 

= t 1 + t 2 + t 3 + u 

where 

t 1 = {n?Tf-T h f,E h g-g) 

t 2 = (n h v Tf-T h f,g) 

t 3 = {T h (f-nZf),E h g-g) 

U=(T h (f-HZf),g), 

and proceed to estimate the ij's, beginning with t\. Let r 1 = min(s M , k + 1). We use 
Theorem O and (PJ) to get 

tx = (/7?V - u[, (E h - E)g) 

< C T h ri (\u f \ H r 1{n) + \q f \ H n { n))h ri \\g\\ 

<cu 2ri il/u hi 

by <m>- 
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Next, for t 2) w e use a previously known duality identity [TO] Lemma 4.1] by which 

t 2 = (i7*V - u{, g) = (c - <?/) , 77^ 9 - q 9 ) + (q f - 1%?', W - Vw h )^ h 

for any Wh G Wh- The first term on the right hand side can be bounded by C T h 2ri \\f\\ \\g\ 
using (iTTj) and Theorem 12.11 For the second term, we use the standard result 

inf ||Vw 9 - Vw h \\ < Ch ro \u 9 \ H r 0+ i (Q) 

w h &W h v ' 

with r = min(s /J , fc), to get 

* 2 < a/> 2ri ii/ii ibn + c r ^ 1 (k / l^(^) + kV^)) ^°l^l^o +1( c) 

< L7 T (/^ 2ri + /i ri+ro ^ 
by the regularity assumption ( 12"U|) . 



Now consider £3. Since T is a bounded operator in L (fi), by Theorem 13.11 we have 

||T^|| < C T |MI VueL 2 (ft). (30) 

Hence, 

£ 3 = (T h (f - Off), (E h - E)g) < C T \\f - n?f\\ Ch n \\g\\ 

by (UgD and (EOD- To control / - ilf / let us first note that = fif since / G R(E). 
Hence absorbing the /i-dependence into C T , we can write 

||/ " Kf\\ < C T h^(\uf\ H r l{Q) + \q%r l(n) ). (31) 

Using the regularity assumption (1211 again, we conclude that 

t 3 <^ 2ri |l/ll \\9\l 



Finally, consider t 4 . By Lemma [2 .2\ is self-adjoint, so 

= (f-Kfy h -<, k -i) (32) 

where G L 2 (Q) is the solution of the HDG method using polynomials of degree 
k — 1 in place of for all > 1, while when k = 0, we set w^-i = 0- We may introduce 
Mfc_i in (E2D because of ( ITOb]) . By Theorem O and (120]) . 



IK-<fe-ill < \K~ U \\ + < C T h ro \\g\\. 

Hence, using also ( 13T|) . we find that 

t 4 <L7 T ^ 0+ri ||/|| 

Combining the estimates for all ti to bound the right hand side of (129]) . we find from (128]) 
that 



\\(T-T h )f\\ <C T h r ^\ 
Using this in (|27j) the proof is finished. □ 
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5. The hybridized eigenproblem 

The main advantage the HDG method possesses over other DG methods, in source 
problems, is that one can eliminate (all) the interior variables {u h , q[) to obtain a single 
equation for the Lagrange multiplier r] h . Since rf h is defined on mesh faces of dimen- 
sion n — 1, the reduced system for r) h can be significantly smaller in size for high degrees k. 
It is natural to ask if this reduction in system size can be carried over to the eigenprob- 
lem. In this section we study this issue. The analysis here is modeled after [9]. 

First, let us review the above mentioned elimination result for the source problem. 
Define local solution operators Q : Mh H- Vh, IX : Mh h-> Wh, Qw '■ L 2 (Q) H- Vh, 
Uw : L 2 {Vt) i — y Wh, element by element, by the following two systems: For any K in 

(cQn,r) K - (IX/i, W-r) K = -(/i,r ■ n) dK for all f G V h , (33a) 

(w, V • QfJ^K + (r(ll/i — fi),w)gK = for all w G Wh, (33b) 

and 

(c Q w f, r) K - V • f) K = for all f G V h , (34a) 

(w, V • £W/k + (tIW, ™W = (/, ™k for all w G W h . (34b) 
The properties of these operators are amply discussed in [HI E] and will not be repeated. 
Let a h (r],fi) = (cQrj, Qfi),? h + (r(Urj - rj), (It/x - fi))ap h and = (f,U[i)r h . 

Theorem 5.1 (The reduced source problem - see [HI Theorem 2.1]). The functions q[ G 
Vh, u h G Wh, and r] h G satisfy (EJ) z/ and onZy z/ 77^ is £/ie unique function in Mh 
satisfying 

a h (riL ^) = h h{v) for all fi G M h , (35) 

q[ = Q V f h + Q w f and u{ = U V f h + V^f. (36) 

Our aim is to obtain an analogue of Theorem 15.11 for the eigenproblem. To this end, 
we will need the next lemma. Let 

\K\ 

\W\- 

We assume throughout this section there is a mesh-independent > such that 

(dhr) 2 h K < C* and cf r < C*. 
This certainly holds for the most commonly used choice of r = 1. 

Lemma 5.2. The local solution operators satisfy the following bounds: There are positive 
constants Co and C\ ( independent of hx ) such that 

IllWlk < C (dl) 2 h 2 K \\f\\ K V / G L 2 {K) (37) 

lllC^llx < C lC lMh, K y^eM h . (38) 

Moreover, whenever < k < l/(Co(d^ T ) 2 h 2 K ) ; the operator I — nllw is invertible and 

w(i - /d^rvik < — „) JK vl2 11/iu v / e L2 w- ( 39 ) 



J-K 1 1 Z^maxr N-l/2 i , / * t \l/2 ||,,||2 ||,,||2 



1 - KC (dl) 2 h 
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Proof. The estimates fl37|) and f]38]) are proved in [7]. To prove fl39|) . we use f )37|) . by 
which 

< Tll/lk- 

where 7 = K,(d^ T ) 2 Coh'j <: . By the given assumption 7 < 1, so the L 2 (f2)-operator norm 
of kUw is less than one. Consequently, / — kMw is invertible and the norm of the inverse 
is less than (1 — 7) -1 . □ 

Theorem 5.3. Let < C3 < (CoC*) -1 and let \. be any positive number less than 
Cz/h. Then I — XhXkv is invertible, and moreover, Xh satisfies 

a h (r) h , /i) = X h ((/ - X h U w )~ 1 '\lr) h , U/i) Sh V/i e M h (40) 

with some nontrivial rjh in Mh, if and only if the number Xh and the functions 

rj h , u h — (I — XhlkvV^rjh, and q h = Qrj h + X h %-u h (41) 

together solve the HDG eigenproblem (1141) . 

Proof. The argument proceeds as in [9], so we will be brief. Setting / = XhUh in and 
applying Theorem 15. 1] (1361) . we get Uh = Ur/h + Uw(XhUh). Hence, 

u h = (/ - XM^Urj,,, (42) 

where the inverse exists by Lemma [5T21 whenever X h < l/(Co(d^.) 2 /i^). Note that the 
later inequality is satisfied whenever Xh < C 3 /h because XhCod^h^ < XhC C*h K < 
XhC C*h < C 3 C C* < 1. Now, using (142]) in the right hand side of (135]) and proceeding 
as in [9], the proof is finished. □ 

Theorem 15.31 is the analogue of Theorem 15.11 for the eigenproblem. It shows, roughly 
speaking, that the reduced form of the eigenproblem, namely ( 1401) . does not lose eigen- 
values up to 0(l/h). In particular, the physically important lower range of the spectrum 
is preserved. The difficulty with ( 140]) is that in spite of being a smaller system ( 1T41 . it 
is a nonlinear eigenvalue problem. Fortunately, the result we present next shows that 
good initial guesses for the nonlinear eigenproblem P0]) can be calculated by solving a 
more standard symmetric generalized eigenproblem. This standard eigenproblem is to 
find an eigenvalue > and a corresponding eigenfunction rjh ^ in Mh satisfying 

a h (fj h ,ii) = X h (Ufj h ,UiS)*r h V/i <G M h . (43) 

and algorithms for solving P3j) are well developed. Although Xh may not equal A^, it is a 
good approximation to Xh (see Theorem I5.4[) and hence can be used as an initial iterate 
for iterative algorithms (such as those studied in [9]) for the nonlinear eigenproblem. 

Theorem 5.4. LetO < C3 < (C0C*) -1 and consider any eigenvalue Xh of (114]) satisfying 
Xh < Cs/h. Then there are two constants, h > (depending on Xh) and C > 
(independent of Xh), such that for all h < h , there is an eigenvalue Xh of (143]) satisfying 



(44) 
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Proof. The idea, as in [9], is to compare two operators which have Xh and Xh as eigen- 
values. The operators are : M h i— > M h and S^: i— > M h defined by 

^(Sf V,7) = ((/ - «IW) _1 ^,U 7 )^ V 7 G M h , 

a fc (5 h //, 7) = (U/x, 1X7)^ V 7 G Mh. 

We will use S h with k = A^, noting that the inverse of I — Xhliw appearing in the defini- 
tion of S^^ exists since A^ < C3//1 (due to Lemma I5T21 - see also proof of Theorem 15. 3p . 
Both operators are self-adjoint in the %(-, -)-inner product: The self-adjointness of Sh 
is obvious. It is also easy to see that (f,U w g)K = (cQ w g, Q w f)i< + (rUwf, Uw-g)dK due 
to f[3"3"j) . so Uw in self-adjoint in the L 2 (n)-inner product. Therefore, S^ is self-adjoint 
in the ah(-, -)-inner product. 

Now, (1401) and (1431) imply, respectively, that 

A, 1 G a(S { h Xh) ) and A^ 1 G a(^), 



A h A h 



hence, by Weyl's theorem [19] on eigenvalues of self-adjoint operators, we conclude that 

/nc(^) 5 11 _ a fe((^l A,l) ~ S h )l,fi) ... 

<\\S h '-S h \\ a = SUp r-^ \T7o- ( 4 5) 



In the remainder of the proof, we bound the right hand side. 

By Lemma [5.21 choosing ho appropriately small, we have, for all h < ho, 

a h ((S h Xh) - S h )j,fi) = (AVM7 - \ h 'U w )~ 1 U'y,U(i) 

^ 1 Xk ?rr K h H^ll H^il ^ and » 

1 — AhLsoU*llK 

< CX h h\\Ui\\ \\Ufi\\ 

< CA^lMUHHU by (EHD. 

where ||/i|| 2 = Hz-tllh We now use the Poincare-type estimate [7J Theorem 3.4], 
\\nWh < C°/i(a*)A*) to conclude that 

a,((S[ Ah) - 5,)7,/i) < CA h /i a h ( 7 , 7 ) 1 / 2 ah(Ai,/i) 1/2 . 
Returning to (145]) and using this estimate, the theorem is proved. □ 



6. Numerical experiments 

In this section, we present numerical results to illustrate the theoretical results of 
the previous sections. We consider model eigenproblems on a square and an L-shaped 
domain and compute the spectral approximations using the HDG discretization. The 
model problems are the same as those considered in [9] so as to facilitate comparison with 
the HRT mixed method. We also present a local postprocessing technique that enhances 
the eigenfunction and eigenvalue accuracy beyond the convergence orders predicted by 
the theory. We begin by describing this postprocessing in the next subsection and present 
the numerical results in the later subsections. 



14 



J. GOPALAKRISHNAN, F. LI, N.-C. NGUYEN, AND J. PERAIRE 



6.1. Local postprocessing. To postprocess the eigenf unction, we are motivated by 
the theoretical results that show that the approximate gradient converges at the 
same order as the approximate eigenfunction Uh- Accordingly, following [18], we define 
(element by element) the postprocessed eigenfunction u* h G P k+ i(K) by 

(V<, Vw) K = -(cq h , Vw) K , Ww G P k +x(K), (46a) 
«, l) K = (u h , 1) K , (46b) 

for all elements K G 2?h- A convergence theory for this postprocessing (that predicts 
that u* h converges at the rate 0(h k+2 ) for k > 1) is available for solutions of the source 
problem [TS] and for the HRT mixed eigenproblem [S]. Next, we define a postprocessed 
eigenflux q£ as the unique element of [Pk(K)] n + xP k (K) satisfying 

((€ ~ Qh) ■ n, v)f = 0, V/i G P k (F), VF C dK, 

(«& " Vh,v) K = 0, We [P k -i{K)] d if k > 1, 



for all elements K G i^. Note that is i/(div, f2)-conforming. 

Using this postprocessed eigenfunction and eigenflux, we are now motivated by the 
Rayleigh quotient to define the following expression for computing an approximate eigen- 
value: 

y = (aV<, + (g • n, uj) aTh 

As the numerical results below indicate, this postprocessed eigenvalue X* h can be a su- 
perior approximation than A^. 

6.2. The square domain problem. We consider the domain Q = (0,ir) x (0,7r). 
In this case, the exact eigenvalues and eigenfunctions are given by X mn = m 2 + n 2 
and u mn = sm(mx) sin(nx), respectively, for m,n G N+. Clearly, the eigenfunctions 
are infinitely smooth, so the convergence rates should be limited only by the degrees 
of the approximating polynomials. We obtain an initial mesh by subdividing Q into 
a uniform grid of 4 x 4 squares and splitting each square into two triangles by its 
positively sloped diagonal. Successively finer meshes are obtained by subdividing each 
triangle into four smaller triangles. The mesh at "level £" is obtained from the initial 
mesh by i refinements. We compute the solution of the HDG eigenproblem in each of 
these meshes. The results obtained are collected below. 

In Table [U the error and order of convergence of the approximate eigenvalues for 
r = 1 is presented. We see that the approximate eigenvalues A^ converge to the exact 
values at the rate of 0{h 2k+1 ). This is in good agreement with the theoretical result of 
Theorem 14.11 To compare with the HRT result, see [H Table 1], which shows that the 
HRT eigenvalues converge faster at the rate 0(h 2k+2 ). 

Curiously however, we found that the simple local postprocessing fT4Hl) can give eigen- 
values competitive with the HRT eigenvalues. In Table El we present the error and order 
of convergence of the postprocessed eigenvalues Aj*. There is no change in the conver- 
gence order for k = as both the approximate and postprocessed eigenvalues converge 
linearly. However, observe that when k > 1, the postprocessed eigenvalues converge at 
a faster rate of 0(h 2k+2 ). Presently, we do not have a rigorous proof for this conver- 
gence rate. This seems to be the fastest rate we can expect. Indeed, we also observed, 
in unreported experiments, that when the postprocessing (148]) is applied to the HRT 
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method, no further improvement in the rate of convergence beyond the 0(h 2k+2 )-ra,te. 
To conclude the discussion on the postprocessing, see Table El where the error and order 
of convergence of the approximate and postprocessed eigenfunctions is presented. We 
see that the convergence rate of the approximate eigenfunctions is 0(h k+1 ), while the 
convergence rate of the postprocessed eigenfunctions is 0(h k+2 ) for k > 1 and 0(h k+l ) 
for k = 0. These results illustrate that the local postprocessing is effective for k > 1 as it 
increases the convergence rate of both the eigenvalue and eigenfunction approximations 
by one order. 

In Table HI we display the absolute value of difference between the approximate eigen- 
values A/j and the perturbed eigenvalues A^ of Section 5, namely \\ h — \ h \. We see that 
the numbers A^ can serve as good approximations of A^. The difference |A^ — A^l is seen 
to decrease with h at the rate 0(h) for k = 0. This is in accordance with the estimate 
of Theorem 15.31 However, the same difference decreases at the rate 0(h 2 ) for k > 1, 
which is one order higher than that predicted by the theory. These results show that 
by solving a standard generalized eigenproblem for A^, we obtain very effective initial 
iterates for nonlinear iterative algorithms to solve the condensed nonlinear eigenproblem 
for A^. At the same time, we note that A/j does not converge to the exact eigenvalue A 
at as fast a rate as A^. 

Finally, we examine the performance of the HDG method for different choices of r. In 
particular, we show the error and order of convergence of the approximate eigenvalues 
for r = h in Table [5] and for r = 1/h in Table |6j We observe in both cases that the 
approximate eigenvalues converge at the rate 0(h 2k ). This convergence rate is one order 
less than the convergence rate of the approximate eigenvalues for r = 1. A theoretical 
explanation for this phenomena follows from ( II ip . In each of the case r = h and t = 1/h, 
one of the bounds in (jlip deteriorate by one order in h. When these revised estimates are 
used in the ensuing eigenvalue convergence analysis, we obtain the reduced 0(h 2k )-rate. 



degree 


mesh 


first mode 


second mode 


fourth mode 


sixth mode 


k 


I 


error 


order 


error 


order 


error 


order 


error 


order 







4.49e-2 




1.26e-0 




3.03e-0 




4.58e-0 






1 


3.45e-2 


0.38 


7.11e-l 


0.82 


1.86e-0 


0.70 


3.00e-0 


0.61 





2 


2.05e-2 


0.75 


3.82e-l 


0.90 


1.07e-0 


0.79 


1.78e-0 


0.75 




3 


l.lle-2 


0.89 


1.99e-l 


0.94 


5.84e-l 


0.88 


9.81e-l 


0.86 




4 


5.74e-3 


0.95 


1.01e-l 


0.97 


3.06e-l 


0.93 


5.16e-l 


0.93 







5.97e-3 




1.40e-l 




7.35e-l 




1.77e-0 






1 


8.44e-4 


2.82 


1.58e-2 


3.14 


9.52e-2 


2.95 


2.36e-l 


2.91 


1 


2 


1.10e-4 


2.94 


1.84e-3 


3.11 


l.lle-2 


3.10 


2.65e-2 


3.15 
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1.39e-5 


2.98 


2.20e-4 


3.06 


1.32e-3 


3.07 


3.08e-3 


3.10 
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1.75e-6 


2.99 


2.69e-5 


3.03 


1.61e-4 


3.04 


3.71e-4 


3.06 







1.38e-4 




3.65e-3 




4.15e-2 




1.27e-l 






1 


4.53e-6 


4.93 


1.04e-4 


5.13 


1.26e-3 


5.04 


3.51e-3 


5.18 


2 


2 


1.43e-7 


4.98 


3.12e-6 


5.06 


3.82e-5 


5.05 


1.02e-4 


5.10 




3 


4.50e-9 


4.99 


9.52e-8 


5.03 


1.17e-6 


5.03 


3.06e-6 


5.06 




4 


1.41e-10 


5.00 


2.94e-9 


5.02 


3.60e-8 


5.02 


9.37e-8 


5.03 



Table 1 . Error and order of convergence of the approximate eigenvalues 
Xh for r = 1. 
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0.50 


o 


2 


1.09e-l 


0.90 
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0.70 
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3.40e-l 


0.91 
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4 
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1 
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8.95e-6 


4.08 


4.88e-5 


4.39 


2.56e-4 


4.69 


3.79e-4 


4.98 
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5.42e-7 


4.05 


2.65e-6 


4.20 


1.18e-5 


4.44 


1.57e-5 


4.60 




4 


3.33e-8 


4.02 


1.54e-7 


4.10 


6.17e-7 


4.25 


7.73e-7 


4.34 







5.40e-5 




7.03e-4 




1.06e-2 




2.14e-2 
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8.04e-7 


6.07 


9.39e-6 


6.23 


1.43e-4 


6.20 


2.14e-4 


6.64 


2 
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l.lle-8 


6.18 


1.09e-7 


6.43 


2.01e-6 


6.16 


2.73e-6 


6.29 




3 


1.28e-10 


6.44 


1.32e-9 


6.37 


2.74e-8 


6.20 


3.34e-8 


6.35 




4 


1.56e-12 


6.35 


2.19e-ll 


5.91 


3.3e-10 


6.38 


2.88e-10 


6.86 



Table 2. Error and order of convergence of the postprocessed eigenvalues 
A* for r = 1. 
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£ 


\\ u ~ u h\\.T h 

error order 


error order 


l n ~ u h\\s h 

error order 


lk-<lk 

error order 






1 

2 
3 
4 


4.19e-2 

1.97e-2 1.09 
9.46e-3 1.06 
4.62e-3 1.03 
2.28e-3 1.02 


8.65e-2 — 
2.72e-2 1.67 
1.03e-2 1.40 
4.65e-3 1.15 
2.27e-3 1.04 


5.55e-l — 
1.68e-l 1.73 
6.83e-2 1.30 
3.09e-2 1.14 
1.47e-2 1.07 


6.11e-l — 
1.91e-l 1.68 
7.27e-2 1.40 
3.14e-2 1.21 
1.47e-2 1.09 


1 



1 

2 
3 
4 


5.02e-2 — 
1.26e-2 2.00 
3.14e-3 2.00 
7.85e-4 2.00 
1.96e-4 2.00 


8.84e-3 — 
1.14e-3 2.96 
1.44e-4 2.98 
1.82e-5 2.99 
2.28e-6 2.99 


3.21e-l — 
8.10e-2 1.99 
1.92e-2 2.07 
4.67e-3 2.04 
1.15e-3 2.02 


7.84e-2 — 
9.23e-3 3.09 
1.13e-3 3.03 
1.39e-4 3.02 
1.72e-5 3.01 


2 



1 

2 
3 
4 


5.28e-3 — 
6.70e-4 2.98 
8.41e-5 2.99 
1.05e-5 3.00 
1.32e-6 3.00 


8.36e-4 — 
5.29e-5 3.98 
3.32e-6 4.00 
2.07e-7 4.00 
1.30e-8 4.00 


6.98e-2 — 
8.10e-3 3.11 
9.78e-4 3.05 
1.20e-4 3.02 
1.49e-5 3.01 


1.31e-2 — 
7.26e-4 4.17 
4.37e-5 4.06 
2.70e-6 4.02 
1.68e-7 4.00 



TABLE 3. Error and order of convergence of the approximate (uh) and 
postprocessed eigenfunctions for r = 1. 



6.3. The L-shaped domain problem. To study the limitations imposed by singu- 
larities of eigenfunctions, we consider the L-shaped domain Q = Q \Qi, where Qq = 
(0, 2) x (0, 2) and Qx = (1, 2) x (1, 2) are the square domains. This domain has both sin- 
gular and smooth eigenfunctions, so offers an interesting example to study the changes 
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4 


2.03e-4 


2.02 


1.30e-3 


2.05 


3.44e-3 


2.10 


5.48e-3 


2.12 







3.23e-2 




2.09e-l 




5.86e-l 




9.78e-l 






1 


8.11e-3 


1.99 


5.10e-2 


2.03 


1.32e-l 


2.15 


2.09e-l 


2.23 


2 


2 


2.04e-3 


1.99 


1.28e-2 


2.00 


3.28e-2 


2.02 


5.13e-2 


2.03 




3 


5.12e-4 


1.99 


3.20e-3 


2.00 


8.20e-3 


2.00 


1.28e-2 


2.00 




4 


1.28e-4 


2.00 


8.02e-4 


2.00 


2.05e-3 


2.00 


3.21e-3 


2.00 



Table 4. The absolute differences |A^ — A/J and their rate of convergence 
to for r = 1. 



degree 


mesh 


first mode 


second mode 


fourth mode 


sixth mode 


k 


£ 


error 


order 


error 


order 


error 


order 


error 


order 







1.42e-l 




1.72e-0 




3.76e-0 




5.45e-0 






1 


3.74e-l 


-1.39 


2.01e-0 


-0.22 


4.15e-0 


-0.14 


5.82e-0 


-0.09 





2 


4.33e-l 


-0.21 


2.09e-0 


-0.05 


4.27e-0 


-0.04 


5.91e-0 


-0.02 




3 


4.48e-l 


-0.05 


2.11e-0 


-0.01 


4.30e-0 


-0.01 


5.93e-0 


-0.01 




4 


4.52e-l 


-0.01 


2.11e-0 


0.00 


4.31e-0 


0.00 


5.94e-0 


-0.00 







1.29e-2 




2.45e-l 




1.31e-0 




2.47e-0 






1 


4.44e-3 


1.54 


9.15e-2 


1.42 


5.57e-l 


1.23 


1.24e-0 


0.99 


1 


2 


1.38e-3 


1.69 


2.55e-2 


1.84 


1.76e-l 


1.66 


4.64e-l 


1.42 




3 


3.62e-4 


1.93 


6.56e-3 


1.96 


4.73e-2 


1.90 


1.37e-l 


1.76 




4 


9.15e-5 


1.98 


1.65e-3 


1.99 


1.20e-2 


1.97 


3.60e-2 


1.93 







1.81e-4 




6.85e-3 




7.60e-2 




2.19e-l 






1 


1.57e-5 


3.53 


6.21e-4 


3.46 


7.18e-3 


3.40 


1.89e-2 


3.54 


2 


2 


1.25e-6 


3.65 


4.24e-5 


3.87 


4.93e-4 


3.87 


1.26e-3 


3.90 




3 


8.23e-8 


3.92 


2.70e-6 


3.97 


3.15e-5 


3.97 


8.01e-5 


3.98 




4 


5.21e-9 


3.98 


1.70e-7 


3.99 


1.98e-6 


3.99 


5.02e-6 


3.99 



Table 5 . Error and order of convergence of the approximate eigenvalues 
Xh for r — h. 



in convergence rates due to singularities. As before, we consider triangular meshes that 
are successive uniform refinements of an initial uniform mesh. The initial mesh is ob- 
tained as in the previous example using a 4 x 4 uniform grid of Qq, except we now omit 
all triangles in f^. 

Since Q has a reentrant corner at the point (1, 1), some eigenfunctions are singular. 
Specifically, we may only expect Theorem 4.1 to hold with s\ = | — e for an arbitrarily 
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flprrrpp 


mesh 


first mode 


second mode 


fourth mode 


sixth mode 


k 


i 


error 


order 


error 


order 


error 


order 


error 


order 




o 


2.53e-l 




7.42e-l 




2.19e-0 




3.57e-0 






1 


5.73e-l 


-1.18 


7.76e-l 


-0.07 


8.88e-l 


1.30 


6.80e-l 


2.39 


o 


2 


6.68e-l 


-0.22 


1.33e-0 


-0.78 


2.18e-0 


-1.29 


2.75e-0 


-2.01 




3 


6.92e-l 


-0.05 


1.49e-0 


-0.16 


2.56e-0 


-0.23 


3.39e-0 


-0.30 




4 


6.99e-l 


-0.01 


1.53e-0 


-0.04 


2.66e-0 


-0.05 


3.56e-0 


-0.07 







1.46e-2 




5.70e-2 




3.88e-l 




2.92e-l 






1 


5.84e-3 


1.32 


1.95e-2 


1.55 


5.78e-2 


2.75 


9.51e-2 


1.62 


2 


2 


1.59e-3 


1.88 


6.75e-3 


1.53 


2.33e-2 


1.31 


3.05e-2 


1.64 




3 


4.06e-4 


1.97 


1.80e-3 


1.91 


6.36e-3 


1.87 


8.70e-3 


1.81 




4 


1.02e-4 


1.99 


4.57e-4 


1.98 


1.62e-3 


1.97 


2.24e-3 


1.96 







2.78e-4 




1.12e-3 




1.76e-2 




2.47e-2 






1 


2.55e-5 


3.45 


1.82e-4 


2.62 


l.lle-3 


3.99 


2.10e-3 


3.57 


2 


2 


1.72e-6 


3.89 


1.51e-5 


3.59 


1.02e-4 


3.44 


1.42e-4 


3.87 




3 


1.10e-7 


3.97 


1.00e-6 


3.91 


6.89e-6 


3.89 


9.97e-6 


3.83 




4 


6.88e-9 


3.99 


6.34e-8 


3.98 


4.39e-7 


3.97 


6.40e-7 


3.96 



Table 6 . Error and order of convergence of the approximate eigenvalues 
Ah for r = 1/h. 



small e > for singular eigenfunctions. For this L-shaped domain the first eigenmode is 
singular and the corresponding eigenvalue is calculated in p[| as Ai = 9.63972384464540. 
It is interesting to note that the third eigenmode is smooth and the third eigenvalue is 
known exactly as A3 = 2tt 2 . 

The errors and resulting order of convergence for the approximate and postprocessed 
eigenvalues are reported in Table [7] for the first and third eigenmodes. We observe 
that the convergence rate of the approximate smallest eigenvalue is at most 0{h 4 ^) 
which agrees with the a priori error estimate given by Theorem 4.1. Furthermore, 
the postprocessed smallest eigenvalue also converges at the same order 4/3 for k > 
1. However, for the third eigenmode, the approximate eigenvalue converges at order 
0{h 2k+1 ) and the postprocessed eigenvalue converges at order 0{h 2k+2 ) for k > 1. 
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eigenmode 


first 


third 


degree 
k 


mesh 

£ 


IA - X h \ 
error order 


IA - Ail 
error order 


IA - AJ 
error order 


IA - Ail 
error order 
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2 
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